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ABSTRACT 

Experimental spectra and images of the supernova remnant SN 1006 have been re- 
ported for radio, X-ray and TeV gamma-ray bands. Several comparisons between mod- 
els and observations have been discussed in the literature, showing that the broad- 
band spectrum from the whole remnant as well as a sharpest radial profile of the 
X-ray brightness can be both fitted by adopting a model of SN 1006 which strongly 
depends on the non-linear effects of the accelerated cosmic rays; these models predict 
post-shock magnetic field (MF) strengths of the order of 150 fj,G. Here we present 
a new way to compare models and observations, in order to put constraints on the 
physical parameters and mechanisms governing the remnant. In particular, we show 
that a simple model based on the classic MHD and cosmic rays acceleration theories 
(hereafter the 'classic' model) allows us to investigate the spatially distributed charac- 
teristics of SN 1006 and to put observational constraints on the kinetics and MF. Our 
method includes modelling and comparison of the azimuthal and radial profiles of the 
surface brightness in radio, hard X-rays and TeV 7-rays as well as the azimuthal vari- 
ations of the electron maximum energy. In addition, this simple model also provides 
good fits to the radio-to-gamma-ray spectrum of SN 1006. We find that our best-fit 
model predicts an effective MF strength inside SN 1006 of 32 /iG, in good agreement 
with the 'leptonic' model suggested by the HESS Collaboration (2010). Finally, some 
difficulties in both the classic and the non-linear models are discussed. A number of 
evidences about non-uniformity of MF around SN 1006 are noted. 
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1 INTRODUCTION 

The supernova remnant (SNR) SN 1006 is one of the most 
interesting objects for studies of Galactic cosmic rays. It is 
quite symmet rical with a ra ther simple bilateral morphology 
in radio (e.g. [ Petruk et al.l [2009c'). nonthermal X-rays (e.g 
iMiceli et al.l 120091 ') and TeV 7-rays (|Acero et al.ll201(J ). Its 
prominent feature is the positional coincidence of the two 
bright nonthermal limbs in all these ban ds, including TeV 7 - 
rays as demonstrated by recent results of lAcero et al.ll|201o[) . 

Current investigations of SNRs with TeV 7-ray emis- 
sion demonstrate an ambiguity in the explanation of the 
nature of TeV 7-rays. Namely, the broad-band (radio-to-7- 
rays) spectrum of these SNRs can be fitted by assuming the 
TeV radiation eith er as leptonic or as had r onic in origin (e.g. 
RX J 1713.7-3946: lAharonian et al.ll2006l : iBerezhko fc Volkl 
I2OQ6I ). 

The question of the origin of TeV 7-rays is closely re- 
lated to the problem of the presence and the role of non- 



linear effects of cosmic rays acceleration by the forward 
shock. One of the key parameter distinguishing between 
these two possibilities is the strength (and thus the nature) 
of the post-shock magnetic field. The classical picture con- 
siders only the compression of the typical interstellar mag- 
netic field (ISMF) _Bo ~ 3 fiG to downstream values of the 
order of tens nG. Models including non-linear acceleration 
(NLA) predict that the ISMF is first amplified upstream due 
to the back reaction of accelerated protons to Bo ~ 30 /iG 
and then compressed above hundred /iG. In the former case 
the inverse-Compton (IC) emission of electrons would be re- 
sponsible for most of the TeV 7-rays, in the latter case the 
proton-origin TeV 7-ray radiation is expected to be domi- 
nant. 

The spectrum of SN 1006 may be explained in these two 
scenarios. One limiting possibility (we call it 'extreme NLA 
model'), namely the case of ISMF amplified and compressed 
to Ba ~ 150 /iG is considered in details by IBerezhko et akl 
(I2OO9I ). The model successfully fits the broadband nonther- 
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mal spectrum from SN 1006 and the sharpest radial profile of 
the X-ray brightness. TeV 7-rays are shown to be produced 
in both the inverse-Compton mechanism and the pion-decay 
one, the latter is dominant. 

Here, we present a new method to compare models and 
observations. In particular, we investigate the origin of the 
patterns of nonthermal images in radio. X-rays and 7-rays. 
At present time, this can be done only by using the classic 
MHD and particle acceleration theories. Therefore, the ques- 
tions behind the present paper are: may a classical model 
explain the radio-to- TeV-7-ray observations of SN 1006 and 
can one put observational constraints on some properties of 
the particle kinetics and/or on the MF? 

In this work, we introduce a "classic" model describ- 
ing SN 1006 and compare the spatial distribution of surface 
brightness derived from the model with those from obser- 
vations in different wavelength bands. The comparison will 
allow us to put some constraints on the parameters of the 
model, thus deriving some hints on the physical mechanisms 
governing the cosmic rays acceleration in SN 1006. In the 
following, the section order is determined by the order of 
parameters determination: the azimuthal and radial profiles 
in the radio band are analysed in Sect. (2] the variation of 
the break frequency in Sect. [31 the broadband spectrum of 
SN 1006 is calculated in Sect. [4] to check the consistency of 
our model and to determine the average MF; the X-ray and 
7-ray brightness are investigated in Sect.[Sl Finally, we draw 
our conclusions in Sect. [H] 



2 CONSTRAINTS FROM RADIO MAPS 

We consider an SNR expanding through a uniform ISM and 
uniform ISMF, in the adiabatic stage of its evolution. The 
Sedov solution is therefore appropriate to describe the hy- 
drodynamics of the system. We consider ideal gas with the 
adiabatic index 7 = 5/3. The MF evolution is treated in 
the classic framework, without non-linear amplification. Its 
strength decreases downstrean j^ far away from the shoch 
front and is modelled following iRevnoldd (| 19981 ): the role 
of the ejecta is considered to be negligible. The classic (un- 
modified) shock creates the energy spectrum of relativistic 
electrons in the form N{E)dE = KE'^ exp{~E/Em^,^)dE 
with the spectral index s, normalization K and maximum 
energy of accelerated electrons i?max. Their dependences on 
obliquity are denoted as Ks{Qo) = -R'sH /k(6o), i?max(6o) = 
£'max||/E(6o) where the symbol '||' marks values at the par- 
allel shock. The downstre am evo l ution of the electron spec- 
trum is modelled as in iRevnoldd l|l998D . 

2.1 Azimuthal profiles 

The symmetrical bright limbs in SN 1006 limit the possible 
orientations of the ISMF in the plane of the sky. Close to the 
shock, the azimuthal distribution of the radio surface bright- 
ness Si is mostly determined by Sr{<p) oc <;(<^) crB(v3)'°^^'''^ 



iPetruk et al.l l2009d ') where q is the injection efficiency (a 
fraction of accelerated electrons), (Jb the compression fac- 
tor for MF (unity for parallel and 4 for perpendicular 
shock), if the azimuthal angle. If the injection is isotropic, 
?(6o) = const where Oo is the obliquity angle), then the 
bright radio limbs correspond to projection of the equato- 
rial belt with NW-SE orientation of the ISMF (BarMF or 
barrel- like model). If, on the other hand, the injection prefers 
quasi-parallel shocks, the bright limbs of SN 1006 are two 
polar caps and the ISMF should be oriented in the NE-SW 
direction (CapMF model). 

The ISMF creates an aspect an gle (j>o with the line of 
sight. Recentlv. iPetruk et al.l (|2009d ) have shown that the 
comparison of the experimental azimuthal profiles of the ra- 
dio brightness with those derived from theoretically synthe- 
sized radio images can be a powerful tool to determinethe 
model of electron injection, the MF orientation in the plane 
of the sky and the aspect angle (j>o- In particular, these au- 
thors have shown that, under the assumptions of uniform 
ISMF/ISM, the injection is be isotropic, the 3D morphology 
of the remnant is BarMF and the aspect angle is = 70° 
(these resul ts were confirmed rec ently by detailed MHD cal- 
culations of lSchneiter et ai1l201Gl ). For the sake of generality, 
we explore in the next subsection also the CapMF model 
with the same value of (j}o. 



2.2 Radial profiles 

The post-shock value (denoted hereafter by the index 's') of 
the spectrum normalization Ka is proportional to the injec- 
tion efficiency The injection efficiency may vary with the 
shock strength (velocity). We assume that Kb oc V~'' where 
V is the shock velocity, and & is a parameter. 

In Appendix I All we show that the surface brightness 
distribution of a Sedov SNR in the radio band is 

Sr = const Sr{p,(p;(t>o,b) !/"<"-'>''^A',||B^''+'>''^7? (1) 

where p — p/R, and p is the coordinate along the radius of 
the remnant R. Si accounts for the evolution of the electron 
energy spectrum and MF inside the SNR. For fixed (fi, Si{p) 
is an universal profile which, for a given dependence of ?(0o), 

feet angle and index s, depends only on the parameter 

We use the experimental radio image of SN 1006 pre- 
sented in IPetruk et al.l (l2009d ) to determine the parameter 
b. We extract the radial profiles of radio brightness from the 
regions shown in Fig. [1] The profiles are reported in Fig. [2] 
together with the theoretical profiles Si{p; b) calculated nu- 
merically for three values of 6 = —1, 0, 1. Close to the shock 
front, the experimental profiles seems to be between the the- 
oretical ones calculated for b = —1 and 6 = 0. 

On the other hand, Fig.[2]shows also that the theoretical 
profiles of radio brightness calculated from a Sedov model 
of SNR expanding through a uniform ISM/ISMF do not fit 
the experimental data to larger extent, namely for p < 0.94_R 



^ At the shock front the perpendicular component of the MF is 
enhanced by a factor 4; then the parallel and perpendicular com- 
ponents evolve independently downstream of the shock, following 
the MF flux con servation and flux freezing conditions respectively 
dRevnolddligoi'l . 



^ Eq. UJ shows also that the universal azimuthal profile of the 
radio brightness Si{ip) depends only on the aspect angle ipo be- 
cause b is assumed to be independent of obliquity. This property 
allowed us to determine ipo from the radio map (see Sect. 12.111 
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Figure 1. NE and SW limbs of SN 1006 in radio at A ~ 20 cm (top panels) and X-rays with energy 2-4.5 keV (bottom panels) 
(jPetruk et al.ll2"009d : lMiceli et al.ll2009l V The maximum value of brightness is 100 times the minimum one, in the radio and X-ray images. 
Radio image is smoothed with Gaussian with a = 0.4' to lower fluctuations. Color straight lines mark the regions used for extraction of 
the radial profiles of brightness; length of regions shown is from 0.8R to I.IR. Green lines represent X-ray contours, linearly spaced. 



(see the inset in Fig. [2]). This result may be explained if ei- 
ther the ISMF or the ISM in the neighbourhood of SN 1006 
is not uniform. In fact, the radio brightness is higher where 
either the ISMF strength Bo or the ISM density rio is larger: 

a (s-|-l)/2 r,(s+l)/2 , J- i 

bi cc ^CTjj oc UoBo . As a consequence, a gradient 
of ISMF/ISM can cause vari ous asymetries in the brightness 
distribution of the remnant (jOrlando et al.|[2007l ): e.g. if the 
gradient has a component along the line of sight, it may in- 
crease the brightness insid e the project ion depending on its 
orientation and strength (|Petrukll200ll ). The radio profiles 
from the SW limb support such a scenario: they monoton- 
ically increase from the shock to ~ 0.85-R (Fig. [T| while 
the maximum of the radio brightness in Sedov SNR should 
be located around ~ 0.977?. We investigate the effects of 
a nonuniform ISMF on the remnant morphology in a com- 
panion paper (Bocchino et al. 2010, in preparation), where 
we consider an MHD model of SNR expanding through a 
nonuniform ISMF and compare the synthetic images in the 
radio band with observations of SN 1006. In the present 
study, we adopt 6 = 0. 



3 CONSTRAINTS FROM OBLIQUITY 
DEPENDENCE OF THE MAXIMUM 
ENERGY 

In this section, we aim at deriving some constraints on 
the modeling of SN 1006 from the obliquity dependence 
of the maximum energy deduced from the observations. 
iMiceh et al.) Hooi) considered a set of 30 regions covering 



T 

- b=--\ 




p, a.u. 

Figure 2. Radial profiles of the radio brightness in NE limb of 
SN 1006. Experimental profiles (from regions 1-5, Fig. [TJ are in 
color. Theoretical profiles are in black, for b = —I (dot), b = 
(solid), 6 = 1 (dashed). They are calculated for s = 2, ipo = 70°, 
isotropic injection and ip = 70° (observational profiles are taken 
for ip = 65° — 75°). Internal plot is a zoom-out of the main figure. 



the entire rim of the shell of SN 1006. The spectral fitting of 
the X-ray emission extracted from these regions allowed to 
derive the azimuthal variation of i/brcak, a parameter in the 
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azimuth, degrees from SE 



Figure 3. Variation of -Emax (shown with vertical one-sigma er- 
rors) over the forward s hock in SN 1006, obtained from experi- 
mental data on inbreak (Miccli ct al. 2009) and Eq. l(2]|, for two 
models of ISMF: BarMF (blue crosses), CapMF (red crosses). 
Aspect angle (f>o = 70°. Dashed red line: the loss- limited model 
with Tj = 1, CapMF. Solid blue line: the time-limited model with 
r) = 1.5, BarMF. Dashed blue line: the loss-limited model with 
7) = 7.6, BarMF. 
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Figure 4. Obliquity dependence of the function /e(9o), deter- 
mining the azimuthal dependence of i?max. The blue solid line 
and red solid line correspond to the time-limited and loss-limited 
model of lRevnoT ds (1998). Dotted and dashed lines correspond to 
the same models but using different values for the gyrofactor rj. 



srcut model of XSPEC which is related to the maximum 
energy of electrons as 



-^max — C 



-1/2 1/2 



-1/2 



(2) 



where Bs is the strength of the post-shock MF, ci = 6.26 x 
10^* cgs. We use the above relation together with the exper- 
imental data on Uhreak to determine the azimuthal variation 
of the electron maximum energy -Emax(v)- 

The dependence of Smax on the obliquity angle Oo can 
be represented as Smax(6o) = -Ejnax|i/E(6o) where /e(6o) 
is a smooth function of 0o. The azimuthal profiles of -Emax 
determined with Eq. ([2]) for two different configurations of 
the ISMF (BarMF and CapMF) is shown on Fig. O 

What cause the limitation of -Emax in SN 100 6? In the 
frame work of the classical theory of acceleration, i Reynolds! 
l|l998l ) (to which the reader is referred to for more details) 
developed three different theoretical models for the surface 
variation of -Emax (and, therefore, for its obliquity depen- 
dence). Namely, the maximum energy of accelerated elec- 
trons may be determined: 1) by the electron radiative losse^ 
(in the following loss-limited model), 2) by the limited time 
of acceleration (time-limited model), and 3) by escaping 
of particles from the region of acceleration (escape-limited 
model). The third model results in constant -Emax that con- 
tradicts th e obliqu ity dependence of -Emax derived from ob- 
servations (jMiceli et al. 2009 ; see Fig. [3|. Therefore, we do 
not consider it in the rest of the paper. In the other two mod- 
els, /e(0o) depends basically on the MF compression ratio 
(tb, and on the level of turbulence which is reflected by the 
"gyrofactor" rj > 1, i.e. the ratio between the mean free path 
of p article along the magnetic field and its Larmour radius 
(see lRevnold3ll998h . Figure [4] shows the obliquity depen- 
dence of Emax in the time-limited and loss-limited models 
for different rj. 



^ Unless otherwise stated, we consider radiative losses to be due 
only to the synchrotron emission. IC losses of relativistic elec- 
trons on the cosmic microwave background radiation (CMBR) 
are inefficient to produce prominent changes in our results. 



The obliquity angle is minimum in SN 1006 at azimuth 
0° for BarMF and at azimuth 90° for CapMF. Therefore, 
Emax is expected to increase or decrease with obliquity for 
BarMF and CapMF respectively (Fig. [3]). In the loss- limited 
model of Emax, the function ./e(Oo) increases with increasing 
obliquity for rj > 3 (red lines in Fig. |3| . In the time-limited 
model, the function /e(Oo) increases with obliquity for any 
rj (blue lines in Fig. [4]). In contrast, Fig.|3]shows decrease of 
Emax from ifi = 90° to 180° for CapMF (red crosses). Thus, 
the time-limited model and loss-limited model with 77 > 3 
are not applicable if one considers a polar-caps morphology 
of SN 1006. In the loss-limited case, the fastest decrease with 
obliquity is for rj = 1 but it does not fit the experimental 
profile of -Emax for model CapMF (red dashed line on Fig. [3]). 
To the end, the NE-SW orientation of ISMF (CapMF, polar 
caps) is not able to explain observed azimuthal variation of 
Inbreak, uuder assumptious of uniform ISMF/ISM and clas- 
sic MHD /acceleration. We tried also other aspect angles, 
00 > 50°, either with or without the inclusion of IC radia- 
tive losses. However, the conclusion remains unchanged. 

In the BarMF case (blue crosses), the function /e(6o) 
for SN 1006 may be d etermined by fitt ing the experimental 
data with a model of I Reynolds! (|l998h . The best-fit in the 
time-limited model is reached for rj = 1.5 ± 0.02 {x^/doi — 
12.7, solid blue line on Fig. (31). The best-fit for the loss- 
limited model is for rj — 7.6 ± 0.11 (x^/dof = 25.8, dashed 
blue line) but the shape of the fit does not follow well the 
observed one. 

Thus, the azimuthal variation of inbreak niay be ex- 
plained in the framework of the classic MHD/acceleration 
theories. It limits ISMF orientation to only BarMF configu- 
ration, in agreement with the same conclusion obtained from 
azimu thal fits of the radio surface b rightness (iFletruk et al.l 
!2009j ). The time-limited model of !Revnolds! l|l998h with 
rj = 1.5 is the most appropriate for -Emax(©o); we use 
it in the present paper. In this model, the maximum en- 
ergy of accelerated electrons varies with time very slowly 
(jRevno lds 199^1. We assume therefore that Emax is inde- 
pendent on the shock ve l ocity. Similar conclusions are ob- 
tained bv !Katsuda et al.l l|2010l ): the correlation they found 
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between the X-ray flux and the cut-ofi' frequency is against 
the loss-limited model for -Emax; absence of time variation 
of the synchrotron flux supports assumption about constant 
(in time) maximum energy. 

The solid blue line in Fig. |3] shows that iSmax ~ 
8.5(50/20^0)"^/^ TeV at azimuth ip = 0. Since the aspect 
angle (j>o — 70°, this value of -Emax corresponds therefore to 
the obliquity Qo = 20°. It is smaller at t he parallel shock , 
namely E^^\\ = Q.6UE^^ for tj = 1.5 (|Revnoldslll99i) . 
Therefore, 



E„ 



5.4 



_Bo_ 
20 /xG 



-1/2 



TeV. 



(3) 



The same time-limited model predicts E„ 



= 3.25E„ 



4 CONSTRAINTS FROM TOTAL RADIO, 
X-RAY AND TEV GAMMA-RAY 
SPECTRUM 

In the calculations of the synchrotron spectrum, the self- 
similarity of Sedov solutions allows us to represent the com- 
plex picture of the synchrotron emission from the whole 
SNR (which includes the complicate description of the down- 
stream evolution of the fluid elements, the magnetic field 
and the spectrum of relativistic electrons in the SNR inte- 
rior as well as the full single-electron emissivity convolved 
at each point with the electron spectrum) by a single uni- 
versal constant and a modification factor rjsyn- The former 
is a (reduced) integral of the radio emissivity over the SNR 
volume; the latter reflects the deviation of the X-ray spec- 
trum from the power-law (for more details see Appendix 
IB1|) . r^syn is defined as a ratio of the integral (i.e. from the 
whole SNR) synchrotron flux at a given frequency (e.g. at 
the X-rays) to the power-law extrapolation of the radio flux 
to this frequency; obviously, 77syn < 1. 

In a broad band (from radio to X-rays), the synchrotron 
spectrum of the volume-integrated emission from the whole 
SNR may be represented by f Appendix IB 1|) 



(4) 



where Cr is a constant, and d the distance to SNR. The 
constant ( is different for different models (Appendix IB1|) : 
it is C = 2.68 for s = 2.0 or C = 2.77 for s = 2.1 in our 
reference model of SN 1006 (namely BarMF, with isotropic 
injection and 6 = 0). 

The reduced photon energy is deflned as e = i> = 



J^/!^c(-Emax||,Bo), Uc{E , B) = Cl 

chrotron characteristic frequency: 



9.5 e. 



E„ 



10 TeV / \20^G 



Bo 



sin cj)) E B is the syn- 



(5) 



this 



where e^cv is the photon energy in keV. With Eq. ((3]), 
becomes e = 32.6£kcV- 

The reduced fiducial energy ti = 637 (-Bf t-Emax) 
e y paramet e 
iRevnoldsl fl 

of the importance of radiative losses in modification of the 
high-energy end of the electron spectrum and therefore of 
the X-ray and 7-ray spectra and images: radiative losses are 
essential for ef < 1; if ef > 1, the adiabati c losses are dom - 
inant even for electrons with E ~ -Emax l| Reynolds! 1 1998h . 



one of the key parameter fo r modeling the X-ray and 7- 
ray emission (jRevnoldsllfggi ). The energy et is a measure 
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Figure 5. The modification factors r^syn in X-rays calculated from 
our reference model for two sets of parameters, namely Bo = 
25 /iG and s = 2.0 (dashed blue line), and Bo = 12 and s = 2.1 
(solid green line) . The co rresponding modifi cation factors derived 
from the SUZAKU data llBamba et al.l2008l . black and red crosses 
represent XIS spectra from the front-illuminated CCDs and back- 
illuminated CCD respectively) are also shown for the cases s = 2.0 
(lower) and s = 2.1 (upper). 



With Eq. ((3} and the age t = 1000 yrs, the dimensionless 
fiducial energy at parallel shock is 



Efll = 5.8 



Br 



20 ^iG 



-3/2 



where e{±_ is a^E 



max II 



(6) 



52 times smaller because 



both Bs and £/max are larger at the perpendicular shock. 

The modification factor r]syn{v) shows how the syn- 
chrotron spectrum Fsynii') deviates from the power-law de- 
pendence i/"'""^^'^^. The modification factor is deflned to be 
?7syn = 1 for the radio band; it is effective in the X-ray band 
and rather quickly approaches to unity with v decreasing 

below t^c(E^max, B) . 

In a similar fashion, the spectral distribution of the IC 
emission from the whole SNR is (Appendix IB2P 



(7) 



where rjic is the modification factor and the universal 
constant for IC 7-rays (exact definitions are given in Ap- 
pendix [BH). For our reference model of SN 1006 (BarMF, 
with isotropic injection, and b = 0), (^t ~ 0.81 for s = 2 and 
Ct = 0.79 for s = 2.1. 



4.1 Fit to the radio spectrum 

iMiceU et al.) l|2009l ) measured the radio-to- X-ray photon in- 
dex a = (s — l)/2 of the non-thermal component for each 
of the 30 regions selected to cover the entire rim of the shell 
of SN 1006, finding a ~ 0.5. This value is al most within 
1-a error of the best-fit value a = 0.61q q9 (|Allen et al.l 



|2008| ). obtained for the radio fluxes from SN 1006 at 8 dif- 
ferent radio frequencies (most of the fluxes are from iMilnd 
Il971i) . We consider therefore a = 0.5 as possible choice for 
the spectral index of the synchrotron spectrum. The best-fit 
(x'^/dof = 1.0) for these radio data and fixed a = 0.5 is 
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i^r,obsM = 18.4 GHz) 



-0.5 



Jy- 



(8) 



In addition, we consider also the case a = 0.55 (i.e. s — 2.1), 
a value successluUy used in the broad- band model of th e 
synchrotron and IC spectrum of SN 1006 l|Acero et al.ll201oh . 
The best-fit (x^/dof = 0.56) for the same radio data and 
fixed Q = 0.55 is 



Fr,obs(«^) = 18.1(i^/lGHz)- 



Jy- 



(9) 



4.2 Fit to the X-ray spectrum 

Figure [5] compares the X-ray modification factor derived 
from our reference model with that derived from observa- 
tions for two different sets of parameters {Bo, s). Note that 
the use of the modification factor allows us to avoid uncer- 
tainties in the distance of the remnant, its radius, and the 
density of emitting electrons. The experimental modifica- 
tion factor r;syn,obs is calculated from the SUZAKU X-ray 
spectrum Fx obs of the whole remnant SN 1006 (Fig. 6 in 
iBamba et al.|[200^ 1 as the ratio of the observed X-ray spec- 
trum to the extrapolation of the radio spectrum to X-rays 

/ \ -Fx,obs(^) /ir\\ 
r;Byn,obs(«^) = -t; 7-T- (10) 

This definition together with Eqs. ([5} and ([3} make the mod- 
ification factor 7?syn,obs(i-') esseutially independent on the 
MF. Theoretical 77syn(£; tf y ), for fixed values of s, h, /k(6o) 
and /e(0o) is a function of the reduced fiducial energy ef| 
only. This parameter reflects the efficiency of the radiative 
losses on the evolution of electrons with energies around 
S'max and, therefore, on the shape of the synchrotron X-ray 
spectrum. In SN 1006, it is related to Bo through Eq. 

The strength of the ambient MF Bo = 25 piG together 
with s = 2.0 provide agreement between the X-ray modiflca- 
tion factors derived from our reference model and from the 
observations (Fig. [5] blue line) . A smaller value of the MF 
strength. Bo = 12 ^G, fits the SUZAKU spectrum if s = 2.1 
(Fig. [S] green line). 

The value Bo = 25 /xG is close to t hat found in the 
extreme NLA model (jBerezhko et al.|[2009l) . However, NLA 
model assumes that Bo is compressed by the shock to the 
level B « 150 /xG and such high strength is the same ev- 
erywhere in the SNR volume. In contrast, our model allows 
large values of MF strength only close to the perpendicu- 
lar shock where the MF is highly compressed; as a result, 
the average MF strength in the classic model of SN 1006 is 
smaller than that in the extreme NLA case. 



4.3 Fit to the TeV 7-ray spectrum 

Fig. [6] compares the modification factors rjic derived from 
our reference case and from 7-ray observations. The experi- 
mental modification factor r;-,_obs is calculated from th e TeV 
7-ray spectrum Ky.obs of SN 1006 (|Acero et al.ll2O10l ). It is 
evaluated as a ratio rj^^ohs = ^7,obs/^T of the observed 7-ray 
spectrum to the extrapolation of the Thomson IC spectrum 
to TeV 7-rays. The latter is found from the radio spectrum 



(Ft/f,; 



where the ratio [Ft/F-, 



calculated with Eqs. Q and for 77syn 



^ / thcor 

1. Thus, 



, obs 



- 7, obs 
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(s + l)/2 
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Figure 6. The modification factor r^ic in 7-ra ys derived from our 
reference model (lines) and from the data of lAcero et al.l 1I2OIC1I . 



shaded regions) for the sets of parameters (Bo 
and (So = 25aiG, s = 2.0). 



12AtG, 



2.1) 



The transformation of the observed TeV spectrum Fy_obs to 
the modification factor Jj^.obs depends directly on the mag- 
netic field strength. Note that this is not a new way to esti- 
mate the MF streng th but just a different representation of 
the method used bv lVolk etaP (jioOSi ). 

The TeV 7-ray spectrum is almost restored by the pure 
IC emission in the model with Bo = 12 fiG and s = 2.1 
(Fig. m green line). On the contrary, the model with Bo — 
25 /iG and s — 2.0 (which is supported by the X-ray spec- 
trum as well; see Fig. [5]) does not agree with the TeV spec- 
trum (Fig. [6] blue line) if only the leptonic 7-ray emission is 
considered. Larger values of the MF strength result in the 
requirement of an additional component to fit the TeV spec- 



trum, as it is the case in the NLA model of lBerezhko et al 
( 200S ) or in the mixed or hadronic models of lAcero et al 
( 20101 ). 

Since further in the present paper we consider the 
pure leptonic model for the TeV 7-ray emission, we assume 
Bo = 12 pG. It is worth to note the difference between the 
spectral index stot =2.1 derived in this section for SN 1006 
as a whole and those derived from the regi ons covering 
the remnant edge, and resulting in sioc ~ 2.0 (jMiceli et al.l 
|2Q09). 

In Sect. [5] we shall analyse the azimuthal and radial 
profiles of the surface brightness extracted from regions lo- 
cated quite close to the shock and use s = 2.0, as suggested 
by X-ray observations l|Miceli et al.ll2009l ). We checked the 
role of other s, and our calculations (not reported here) show 
that, for s — 2.1, the profiles of brightness are almost the 
same as those reported here. 

However, our calculations (not reported here) show 
that, for s — 2.1, the profiles of brightness are almost the 
same as those reported here. 
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Figure 7. Azimuthal profile of the X-ray brightness at fixed p 
which corresponds to the maximum of the radial distribution of 
brightness at ip = tt/2. The calculations are done for e = 1.2 keV, 
s = 2 and two values of So, namely 12 /iG (solid green line) 
and 25 iiG (dashed blue line). The exp erimental data ar e taken 
from the hard X-ray image of SN 1006 jMiceh et al.ll2009l) . They 
are derived averaging X-ray brightness along radii within annuli 
centered on the remnant (from 13.8' to 14.8' off the center for 
SE-NE and from 14.4' to 15.2' for SE-SW profile). The profile 
from SE to SW is in red, that from SE to NE is in black. 



5 CONSTRAINTS FROM X-RAY AND 
GAMMA-RAY MAPS 

5.1 X-ray azimuthal profiles 

In Appendix lAll we demonstrate that the distribution of 
the surface brightness of a Sedov SNR due to synchrotron 
emission can be represented as 



Sx = const 5x(/>, p, (p; (t>o, b, e^) eI-,^^K^hBoR. 



(12) 



The universal shape <Sx of the radial (for fixed ip) and 
azimuthal (for fixed p) profiles is determined just by one 
parameter, efy , if s, (po, b as well as the obliquity dependence 
of the injection efficiency and the model for the maximum 
energy of electrons are fixed. In case efy ^ 1 and/or v <^1, 
the role of the radiative losses on the downstream electron 
distribution is negligible and the profiles of the brightness is 
then independent on the fiducial energy: 



Sx = const 5r(p, 0o, &) ? ^" ^^^'^E^^^K^wBoR, 



(13) 



that is the same as Eq. IT} for the radio brightness. On the 
other hand, our calculations show that the role of the evo- 
lution of injection efficiency in time (which is represented 
by 6) is less important for X-rays than for the radio be- 
cause the radiative losses (represented by ety) are dominant 
in determining the downstream distribution of X-ray emit- 
ting electrons. 

In our reference model of SN 1006, even efy is not a 
free parameter, because it is determined through Eq. ((6]) 
by the strength of the MF. Fig. [7| compares theoretical and 
experimental results. Synthesized azimuthal profiles of the 
X-ray brightness agree with the observations though the fit 
is not ideal. Simulations reveal that the strength of the MF 
is not important for the azimuthal variation of the X-ray 
brightness. 



Another possibility to change the azimuthal variation 
of the synchrotron X-ray surface brightness in a model 
is to consider a broader encQ of the electron spectrum 

at the shock, e.g. N{E) oc exp (S/Bmax)"') • The 

value a' — 0.5 ~ 0.6 (|Ellison et al. I I2OOOI . and references 
therein) makes the fit even worse. Really, the azimuthal 
distribution of the brightness is roughly proportional to 

exp ^— (-Bm(v')/£'max(<^))" ^ whcrc Em is the energy of elec- 
trons which give the largest contribution to emission at an 
observed frequency. Fig [7] assumes a' = 1; smaller a' results 
in smaller contrasts between azimuth p = Tr/2 and p — 
that is against of the observations. Larger values of the pa- 
rameter, a' = 1-^2, appeared in the loss-limited model 
(|Zirakashvili fc AharonianI l2007l : ISchure et al] I2OI0I ) might 
increase the contrast but the young age of SN 1006 makes 
the radiative losses ineffective in limitation of iJmax. 

The differences in the synthesized and observed profiles 
might be due to nonuniformity of ISMF and/or ISM: larger 
contrasts of ISMF or ISM density between azimuth 0° and 
90° induced by nonuniformity is obviously able to increase 
contrasts. 



5.2 X-ray radial profiles 

The method for the MF strength estimation from the 
radial profile of the X-ray bright ness is d escribed by 
iBerezhko fc Volkl (|2004 ) (see also Bal iet1l2005l ). The radia- 
tive losses of electrons with energy E is E cc E^B^. These 
losses are less important for electrons emitting in radio but 
they are able to modify effectively the energy spectrum 
of electrons radiating X-rays. As a consequence, the syn- 
chrotron rim in X-rays is thinner than that in radio. The 
idea of the method is that the stronger the magnetic field, 
the larger the radiative energy losses experienced by rela- 
tivistic electrons. This leads therefore to the rapid decrease 
of the spatial distribution of electrons behind the shock and 
to a sharp maximum in the radial X-ray brightness profile. 
From the observational point of view, the thinner the rim in 
X-rays the stronger the magnetic field is expected to be. 

Fig. [H] compares the theoretical profiles 5x(p) with data 
from XMM and Chandra. The simulated distribution with 
Bo = 12 pG (green solid line) fit the XMM data. In our 
model, the MF compression factor is ctb = 3.8 at the az- 
imuth tp = 70°. The post-shock MF is therefore Ba ^ 45 pG 
in both NE and SW limbs. This value could be considered as 
an upper limit for an average MF within the limbs because 
some observed profiles are a bit thicker than the theoretical 
one which is shown by the thick green line. The strength 
Bo = 25 pG (long-dashed blue fine) does not fit XMM the 
radial profiles of X-ray brightness. 

However, th e sharpest Chandra profile (Fig. 4A in 
iLong et al.l l2003h may not be explained by Bs ~ 45 pG. 



^ The broadening of the spectrum is observed in SN 1006 
dRevnoldj [19961 : lEUison et a"ul2000l'l: the observed spectrum is 
fitted with a' = 0.5 — 0.6 jEUison et al.. ,2000, and references 
therein). The reason of the broaden ing should b e related to the 
property of the acceleration process jPetruklbOO^ I and not to the 
inhomogenity o f conditions in di fferent places inside the remnant 
as suggested bv lRevnoldj lll996h . 
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Figure 8. Radial profiles of the X-ray brightness in NE (left) and SW (right) limbs of SN 1006. Experimental XMM-Newton profiles 
(from regions 1-5 and 6-10 resp ectively, Fig.[l]l for the photon energy range 2-^4.5keV are in color. The sharpest Chandra radial profile 
(from Fig. 4A in lLong et al.l |2003. photon energy is 1.2 2keV) is shown by the histogram. Theoretical profiles are shown by the thick 
blue dotted line (for Bo = 25 /xG, smoothed to the Chandra resolution by Gaussian with sigma 0.2"), and by the thick green solid 
line (for 12 /iG, smoothed to the XMM resolution by Gaussian with sigma 2.6"). They are calculated at 1.2 keV photons, for azimuth 
f = 70° in our model of SN 1006, s = 2. Theoretical profile for Bo = 25 ^iG smoothed to the XMM-Newton resolution is shown by the 
long-dashed blue line. 



Our model fits this profile if tiie post-sliock field is ~ 
95 /jG (blue do t ted li ne). The same filament was used by 
iBerezhko e\~^ (j2003l ) to deduce ~ 130 field. Our es- 
timate for the thinnest filament is comparable but lower 
than in the NLA model. The reasons of such discrepancy 
are some differences between our and their models. Namely, 
in our model, MF decreases downstream of the shock while 
the extreme NLA m odel assumes un iform MF. In addition, 
we accept (following lRevnold3 Il998[ ) that accelerated elec- 
trons are confined in the fiuid element while ISerezhko et al.l 
(j2003h include diffusion. 

In general, MF estimated from the radial profile of 
X-ray brightness refiects the local conditions. The quite 
large strength of the downstream magnetic field, B ~ 
130 — 150 /iG, adopted in the extreme NLA model, was 
assu med to be the same e verywhere in the SN 1006 inte- 
rior (IBerezhko et al.ll2009l ). This value is reasonable for a 
thinnest NE filament (Fig. 4A in lLong et al.ll2003l ) as it is 
appa rent from the fitting of the radial profile of X-ray bright- 
ness (jBerezhko et al.l l2003l : iKsenofontov et aL I [2005). How- 
ever, the two close radial profiles are already thicker (Fig. 
4B,C in iLong et all \2QQ± suggesting therefore a smaller 
value of B even around the location of the original sharpest 
filament. 

An effective MF inside SN 1006 (i.e. which may be used 
to represent SNR as a whole) is smaller. Let us consider two 
possible effective field values: the volume average 



BdV 



and the radio-emissivity weighted volume average 



{B)o. = / BPdV 



PdV. 



(14) 



(15) 



(The radio-emissivity weighted volume average is higher 
than the volume average because the emissivity P quickly 
decreases downstream; thus, in calculation of the average. 




Figure 9. Radial profiles of the IC 7-ray brightness at 1 TeV 
in our model of SN 1006, for azimuth 1/3 = (dotted lines) and 
if = ■k/2 (solid lines), for two values of Bo\ (po = 70°, s = 2. 



most of the contribution comes from regions close to the 
shock where the MF is large.) It may easily be shown 
analytically that (due to self-similarity of the Sedov solu- 
tion) both (B)^ and (B)^^ are simply products of Bo and 
some constants. Our model predicts Bo = 12 /iG and yields 
(B)^ = I.O6B0 = 12 /iG and (B)^^ = 2.7Bo = 32fiG. The 
latter is in good agreement with the strength (~ 30 piG) 
found in the leptonic model of lAcero et al.l l|2010l ) and with 
the estima t ion ( again « 30 fiG) derived from Fig. 1 in 
IVolk et all (|2008l ) with the use of the HESS spectrunfl. 



^ Our model deals with three-dimentional distribution of the 
m agnetic field and emitti ng electron s whi le simpler models 
of lAcero et aP l|2010l ^ and IVolk et al.l l|2008l ) consider uniform 
plasma in the uniform magnetic field. 
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5.3 A note on the gamma-ray brightness 

The distribution of the IC brightness in 7-rays with energies 
£ is given by (Appendix I A2[l 

Sic = const 5ic(e,p, V5;0o,6, (;f||,£^max) -ft'sii-R- (16) 

This formula shows which factors affect the shape of the az- 
imuthal and radial profiles and which determine their am- 
phtudes. The 7-ray image of SN 1006 (namely, the shapes 
S of the azimuthal and radial profiles in these bands) - like 
the X-ray map - depend only on the value of Bo, once other 
parameters are fixed. 

SN 1006 is rather faint in TeV 7-rays to allow, at present 
time, to derive azimuthal and radial profiles with quality 
comparable to those obtained in the radio and X-ray bands. 
However, we may check whether our model provides the ob- 
served location of the bright 7-ray limbs. Fig. |5] shows the 
radial profiles of the 7-ray surface brightness in our model 
of SN 1006. The bright limbs are located at the azimuth 
(fi = n/2 for both strengths of MF considered (namely, 12 
and 25/xG), in agreement with the observations. Note that 
such property is not universal; it depends on the parameters 
of the model. For example, if an aspect angle would be 90° 
then the observed location of the TeV 7-ray limbs may be 
possible only for ISMF strength larger than ~ 100 /iG. We 
address this issue in a separate paper. 



6 DISCUSSION AND CONCLUSIONS 

The magnetic field strength in SN 1006 is one of the key pa- 
rameter in the model. Being related to -Emax|| with Eq. ([3]) 
and to the parameter ef|| (which regulates efficiency of the 
radiative losses of relativistic electrons) with Eq. ([§} , it influ- 
ences almost everything in nonthermal spectra and images. 

We consider a 'classic' model of SN 1006, i.e. model 
which is based on classic MHD and acceleration theories. 
Since they are better developed compared to NLA approach, 
they allow us to put observational constraints on the (test- 
particle) kinetics and MF, to compare the azimuthal vari- 
ations of the electron maximum energy and the surface 
brightness in radio, hard X-rays and TeV 7-rays. At the 
present time, such comparison may not be done in the frame 
of the NLA theory. We demonstrate that the 'classic' model 
is in agreement with most of the observational data. 

We try to fix free parameters of the model step-by- 
step, looking for observations which is mostly sensitive to 
some of them (Table [TJ . In addition to the commonly used 
broad-band spectrum, the properties of the nonthermal (ra- 
dio. X-ray and TeV 7-ray) images of SNR as well as spatially 
resolved spectral fits are considered. 

In particular, the morphology and azimuthal profiles of 
the radio brightness may determine the orientation of ISMF. 
Namely, the radio data may be fitted by the model with uni- 
form ISMF which is oriented perpendicular to the Galactic 
plane with an angle 70° to the line of sight (jPetruk et al.l 
l2009d : ISchneiter et"al] I2OI0I ). If so, the injection efficiency 
should be independent of obliquity. The radial distribution 
of the radio brightness depends now only on the way the 
injection efficiency varies with time {K cc V~^). The obser- 
vations however may not definitively fix b. It is somewhere 
between —1 and but accuracy of the data allow also for a 



bit wider range. Spatially resolved X-ray analysis of regions 
around the forward shock demonstrate that distribution of 
Inbreak may be explained by the time-limited model of -Emax; 
this is in agreement with recent results of iKatsuda et al.l 
The maximum energy of electrons at the parallel 
shock is found B^axH = 7(Bo/12 ^G)"^/^ TeV. It is 3.25 
times higher at regions where shock is perpendicular. 

We obtain expressions for the radio. X-ray and 7-ray 
spectra from the whole SNR in a form which clearly show 
which parameter of the model is responsible for the ampli- 
tude of the spectrum and which one for its shape. The mod- 
ification factor of the synchrotron X-ray spectrum - which 
shows the deviation of the spectrum from the power law - 
may well be explained by the classical model with ISMF 
strength Bo = 25^G if s = 2.0 or with Bo = 12 /iG if 
s — 2.1. At the same time, the TeV 7-ray modification fac- 
tor prefers only the pair Bo = 12 ^G, s = 2.1; TeV emission 
is then completely due to IC process. In case Bo = 25 /xG, an 
additional component in the TeV 7-ray spec trum is needed, 
from pion decays, as it is in the model of iBerezhko et all 
(2009). The proton injection in such scenario should increase 
with obliquity in order to fit the observed azimuthal pro- 
files of TeV 7-ray brightness (ISMF is parallel to the limbs). 
Could the electron and proton injections have so different 
dependences on obliquity in the same SNR: isotropic for 
electrons and quasi-perpendicular for protons? 

The extreme NLA approach (|Berezhko et al.l|2009l ) pre- 
dicts Bo — 30 /iG immediately before the forward shock 
and B — 150 fiG everywhere inside the SNR. A number 
of the radial profiles of X-ray brightness obtained from 
XMM image agree with our model if an ambient MF is 
Bo ~ 12fiG. Around the quasi-perpendicular shock, where 
the profiles are extracted from, our model predicts the post- 
shock MF with strength Bs± — 50 /iG. However, in the clas- 
sic model of SN 1006, this is the value immediately post- 
shock; after then it rapidly decreases downstream. There- 
fore, an effective (emissivity weighted average) MF within 
SN 1006 is estimated t o be 3 2 /xG that agrees well wi th es- 
timates of IVolk et ai] (120081 ) and lAcero et al.l l|2010t ). MF 
in the sharpest Chandra profile is fitted in our model with 
Bs = 95 pG; it reflects the local conditions, only within this 
filament. 

We found that the broad-band spectrum from the whole 
SN 1006 is better represented with the electron spectral in- 
dex Stot = 2.1 while local rad io-to-X-ray spectra over the 
SNR shock prefers sioc = 2.0 (|Miceh et al.|]2009l '). It is in- 
teresting, that similar differ ence in spect r al ind ex is found 
in the theoretical study of ISchure et al.l (|2010l ) : the spec- 
trum near the shock is flatter than the overall spectrum. 
The authors attributed this difference to the time evolu- 
tion of -Emax which was lower at previous times. I n contrast, 
the t ime-limited model for the maximum energy l|Revnoldsl 
[1998!) which fits the azimuthal variation of inbreak in SN 1006 
(Sect. [3]) as well as absence of the t ime variation of the 
synchrotron flux (jKatsuda et al.|[2010l ), suggest what i?max 
varies quite slowly in this SNR, at least in the regions close 
to the shock. The issue of different spectral slopes has to be 
considered in the future. As to the purpose of the present 
study, difference between Stot and sioc is negligible for az- 
imuthal and radial profiles of radio. X-ray and IC 7-ray 
brightness. 
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Table 1 . Summary of the observables used for parameter determination and cross-check'' 



Observable 


Parameter 


Value 


h 

radio azumuthal profile 


aspect angle 


J. '7r\o _i_ A 00 
(po = (1) it 4.Z 




injection type 


isotropic 




/^T"! El ^" 1 I~^K^T-r tlTl/H ~i 1\T R m ^~IVT\tl ^^1 ^~imA 
01 iClILcLlilOli Oi lOlVi-i. CL1H_1 OIN-IV lliUl JJliOiUli V 


SF NW barrel like 


raOio racliai pronie 


h \ri K \/~^ 

U Hi -IVg UC V 


— 1 < 6 < 


local broad-band fits of spectra'' 


local index s^q^ over shock 


S2q(. = 2.0 over most of SNR rim 


^reak azimuthal profile'^ 


model of -Bmax 


time-limited"^ 




ratio of the mean free path to Larmour radius 


77 — i.o 




electron maximum energy at parallel shock 


^maxll — '•U(-Oo/li/iCj) ' lev 




electron maximum energy at perpendicular shock 


-£'max_L — ^-^^-^max 1 


radio and hard X-ray spectrum 


MF strength and index stot for the whole SNR 


(Bo = 25 /.iG and stot = 2.0) or 






(_Do = 1^ s>na Stot = ^-Ij 


radio and TeV 7-ray spectrum 


7-ray emission model, MF strength and index stot 


IC with Ba = 12 and stot = 2.1 


X-ray radial profiles 


post-shock MF strength in the limbs 


Bs± ^ 50 AtG 


X-ray azimuthal profile 


MF strength 


OK 




ISMF orientation and aspect angle 


OK 




model of -Bmax 


OK 


7-ray limbs location 


MF strength, aspect angle 


OK 




7-ray emission model 


OK 



the model assumes uniform ISMF/ISM and 7 = 5/3 
^ Petruk et al. (200901 : Schneiter et al. (2010) 
" Miceli et al. (20091 

see also J<atsud^_et_^l^ ^2010) 



Azimuthal profiles of the X-ray and 7-ray brightness in 
our model behave in the same way as in the observations. 

'Classic' model has also few difficulties. 
iRothenfiug et all (|2004 l developed a simple geometrical 
criterion to distinguish between barrel-like and polar-cap 
morphology in SN 1006. They have shown analytically 
that the ratio between the central and the rim brightness 
should be larger than some value in BarMF case (projected 
"barrel" has to provide enough brightness in the internal 
regions). In XMM map, this ratio is smaller. The luck of 
brightness from equatorial belt in the central part is an 
argument against BarMF morphology. Our model, which 
strongly pref er BarMF, do es not agree with the criterion 
of [Rothc nflug et al.l (|2004l '). Nevertheless, the polar-cap 
scenario, which is in agreement with this criterion (and 
is adopted by the NLA model), is unable to explain the 
observed azimuthal profiles of the break frequency inbreak and 
the radio brightness, under assumptions that ISMF/ISM 
are uniform and the amplified/compressed MF increases 
with obliquity. 

Another minor point of the classic model of SN 1006 is 
the rather large ambient MF, Bo ~ 12 /iG, which is difficult 
to expect without MF amplification at the high location of 
SN 1006 above the Galactic plane. 

Our model deals with ideal gas with the adiabatic 
index 7 = 5/3 and cannot explain the small distance 
between the forward shock an d the con tact discontinuity 
ijCassam-Chenai et al.l i2008: M iceh et all [20091. Instead, if 
acceleration is so efficient that relativistic particles affect 
hydrodynamics then the adiabatic index may be smaller 
than ours. The small distance observed may naturally be 
explained by such, more compressible, plasma with the in- 
dex like 7 = 1.1 (|Orlando et al.ll2010l ). It is worth noting 
that in such case of efficient acceleration, there is no need 
for MF amplification: the only shock compression to factor 



a — 21 (as it is for 7 = 1.1) may result in quite large down- 
stream MF even in case of Bo of few /iG. 

The two models, classical and extream NLA, are com- 
pared in Table [21 It is evident that none of them explaine 
the whole set of the SN 1006 properties. A new model of 
SN 1006 has to include either combination of the two ex- 
tremes or inclusion of the ISMF/ISM nonuniformity. 

All results presented here are obtained under as- 
sumption that SN 1006 evolve in the uniform ISMF and 
uniform ISM. It is shown that the scenario of classic 
MHD/acceleration plus uniform ISMF/ISM strongly prefers 
the barell-like morphology of SN 1006. However, we also see 
that nonuniform ISMF/ISM could be an essential element 
in the model of SN 1006. In particular, slanted lobes, the 
inversion of the brightness ratio between NE and SW limbs 
from radio to X-ray band and the higher break frequency in 
NE limb may only be explained by presence of gradient of 
ISMF and/or ISM. We expect that the effect of the nonuni- 
form ISMF might dominate the role of some nonlinear effects 
arising from efficient acceleration of cosmic rays by the for- 
ward shock in SN 1006. We like to address this issue in the 
future. 
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Table 2. Comparison of SN 1006 models 



Properties of SN 1006 



Classic model^ 
(with uniform ISMF) 



Extream NLA model'' 
(with uniform ISMF) 



two-limbs in radio and hard X-ray image 



location of TeV 7-ray limbs 
^othenflug_et_^l^ (2004) criterion 
radio spectrum 
hard X-ray spectrum 
TeV 7-ray spectrum 



radio radial profile 

sharpest X-ray radial profile 

radio azimuthal profile 

hard X-ray azimuthal profile 

TeV 7-ray azimuthal and radial profiles 

i/'break azimuthal profiles 

pre-shock MF strength Bo = 12 /^G 



(eventual) concave shape of the radio spectrum*^ 
very close forward shock and contact discontinuity 
radio 'overbrightness' in the SNR interior 
slanted lobes 

ratio of radio IZ.^ < 1 and X-ray brightness TZ.^ > 1 

'-'break, NE/l'broak.SW > 1 



YES 

SN 1006 is barrel 

ISMF direction: SE-NW 
YES 
NO 

YES as power low 

YES with {B)^^ = 32AtG 

YES 

IC with {B)^^ « 32AtG 
YES 

YES with Bs Ri 95 ^tC 

YES 

YES 

YES 

YES 

NO 

(if ISMF around SN 1006 

is typical ~ 3 ^G) 
NO 
NO 
NO 
NO 
NO 
NO 



YES? 

SN 1006 has polar caps 
ISMF direction: SW-NE 

YES? 

YES 

YES with concave shape 
YES with (B)^ Si 150 AtC 
YES 

IC with (B)^ ss 150 AtG 
and hadronic component 
NO? (uniform B inside SNR) 
YES with Bs » 150 /iC 
? 
7 

7 

NO? 
YES 

as result of amplification 
(if any) 
YES 

YES (with 7 = 1.1) 

NO? 

NO 

NO 

NO 



^ present paper 
fa Berezhko et al. (2009') 
Allen et al. (2008') 
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APPENDIX A: SURFACE BRIGHTNESS OF 
SEDOV SNR 

Surface brightness of a spherical SNR is an integral of vol- 
ume emissivity q along the line of sight 



R 

2 j qdl 





2R 



rvada 



a(p) 



(Al) 



where p is distance from the center of projection, r — r/R, 
a Lagrangian coordinate, — dr/da, 



g= / N{E)p{E,e)dE, 



(A2) 



where E and e are the electron and photon energies, p the 
radiation power of a single electron. In case of Sedov SNR 
in uniform medium th e electron energy distrib ution down- 
stream of the shock is (jPetruk fc Beshlevll2008l ) 



N{E)^KE"£:-/e^p 



E 



-Bmaxlli^adi^rad/E 



(A3) 



the normalization K = Ks\i{t)fK{Qo)K{a), the magnetic 
field B — _Bs|| (t)(7B(0o)-B(a) and the electron maximum en- 
ergy -Emax = £max||/E(6o). 



a(p) 

or, in other form, 

2V3e» (sin 0) A{s) 



^^■2 _ p2 



(A8) 



= n S^VP, <y5 <Po, h) 

race-' 

X (z./ci)-''-^'''^7f.||Br'^/'7?. 



A2 IC emission 

The IC radiation power is 



(A9) 



-E-'l{E,e), 



(AlO) 



where all notations ha ve their common meaning, I is a spe- 
cial integral (see e.g. IPetruklliooi l. The IC brightness is 
therefore 

Sic = 5ic(e,p, v';</'o,&, Cfih-Bmax) KsiiR. (All) 

The function Sic{p,ip) is not so universal as in case of the 
synchrotron emission; it depends on the absolute values of 
the photon energy and the maximum electron energy; we do 
not present it here. 



Al Synchrotron emission 

The synchrotron radiation power is 
3e^ (sin < 



-BE, 



4B 



(A4) 



where all notations have their common meaning. The syn- 
chrotron surface brightness of Sedov SNR is therefore 

2^e^ ( 

6syn = 5 Ssyn{v,p,ip; (po, b, em) E^^^K^^^ BoR.[A5) 

where Ssyn{i',p,i^) is a universal dimensionless function 



a{p) 



X exp 



xaBBfKK 



fadf rad/E 

ffadd 



de 



(A6) 



where e = i?/-Emax|| • It depends on the dimensionless models 
of obliquity variations of K, B and -Emax (i.e. on /k, ctb, /e) 
but is independent of the actual values of -Emax, Ks, Bo and 
R. 

In the limit tf y 2> 1 and/or D <^ 1, Eq. (|A5P transforms 

to 



where 



)A{s 



Sr{p,'p;<po,b) 



(A7) 



APPENDIX B: NONTHERMAL SPECTRUM OF 
SEDOV SNR 



Flux is defined as 
F{u) = (47rd^)"^ / P{v)dV 



(Bl) 



where V is the volume of SNR and P the volume emissivity. 
We assume that the energy spectrum of electrons in the form 

N{E)dE = KE'" exp{-E/Em^^)dE (B2) 

are created at the shock. The volume emissivity is 



P{iy) = / N(E)p(E,u)dE 



(B3) 



where p is the spectral distribution of radiation power of 
'single' electron with energy E. Let us consider adiabatic 
SNR in uniform ISM and uniform ISMF. 

In general, the efficiency of injection may depend on the 
shock obliquity angle Qo- If particles are injected easier at 
quasiparallel shocks then Ks{Oo) is decreasing function of 
Oo with decrement rate d ependent on the le vel of turbu- 
lence, shock strength etc. (jEllison et al.l[l995h . Let us con- 
sider parametric representation Ka = -ft's||/K(0o) with ap- 
proximation /k = exp (— (Go/Ok)'^) where A'ay the normal- 
ization for region immediately after the parallel shock. Ok 
the parameter. Qk = tt/G approximates the classical quasi- 
parallel dependence, ? oc cos^(Oo). In case of the isotropic 
injection, Qk = oo. 



Bl Synchrotron emission 

The radio f lux ||B1[) from Sedov SN R may be written as (for 
details, see lPetruk fc BeshlevllioOTi ') 
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Figure Bl. f for different values of parameters b and Ok- s = 2 



(s+l)/2 n>3 



(B4) 
(B5) 

(B6) 



where 

C= (47r)-U(s)c2M^4^~"/', 
ci = 3e/(47rmcC^), C2 = \/3e^/(m<,c^), 

= (sin((^)'"+^'/^), (C = 3.493 x 10"^" cgs in case s = 2), 
ip the angle between MF and the hne of sight, 



C(6,eK) = / IdesinOfK /dafVai^fcTBB)'"^'^''' (B7) 



o"b(6o) is the compression factor for MF, r and a are Eu- 
lerian and Lagrangian coordinates respectively, ra — dr/da, 
bar represents parameter divided by its post-shock value, 
[ifi, 9) spherical coordinates. Thanks to the self-similarity, 
the constant <^ 'compactifies' th e whole downstream evolu- 
tion of fluid elem ents (j Scdov 1959), magnetic field and rela- 
tivistic electrons (|Revnolds„1998i '). 

In a similar fashion, the X-ray flux is (jPetruk fc BeshlevI 

l2008l ) 



(B8) 



where u = i//i'c(-Emax|| , So), Vc{E,B) oc E^B is the syn- 
chrotron characteristic frequency, C2 = C2 (sin ip) /(47r) a 
constant, ef|| = 637 (-Bg||ii7niax||) is the reduced flducial 
energy. The energy ef is a measure of importance o f radiative 
losses in modification of the electron spectrum l|Revnoldsl 
Il998f ). The function 



Cx(!^;6, 6K,ef||) ^ J '^'^ JddsindfK jdar r-^KaaB 
00 



(B9) 



00 

I' 



X / "C/ exp 



£ad£rad/E 



-Tsyn 




reduced s 



Figure B2. Modification factor r^syn- Calculations are done for 
s = 2, the time-limited model of -Emax with 77 = 1.5, isotropic in- 
jection (thick lines) and quasiparallel injection (thin lines), three 
values of b, efy = 100 (upper panel) and efy = 3.2 (lower panel). 
Experimental modification factor for SN 1006 are shown for com- 
paris on. It is obtained from the SUZAKU spectrum ( Bam ba et al.l 
I2OO8I. Fig. 6) for photon energies > 2 keV, with the use of Eq. mOC . 
MF strength is given by Eq. l[6]l: -Bo = 3 fiG (upper panel) and 
Bo = 30 /iG (lower panel). 



where £ad(ts), frad(i; efy , So) represent adiabatic and radia- 
tive losses of relativistic electrons (Pctruk & Bcshlcy 200^), 
Fsyn the function known in the theory of synchrotron radi- 
ation, e = S/^maxll- 

With i>, the radio fiux (IB4|I may be written in a form 
similar to (IB8I): 



~-(s-l)/2 



ab,QK)KnBoEl-'R^d-\ (BIO) 



Comparison of HB8|) and (|B10|) demonstrates that, for u 
much smaller than X-ray frequencies, (^x transforms to (, 
as expected: 



Cx(i>) = A{s) V 



~-{s-l)/2 



(Bll) 



This transition may also be shown anal ytically from IIB9I). 
in the limit E <: iSmax and E ^ ef-Emax (jPetruk fc BeshlevI 
i200S) . 

Let us introduce the modification factor for the syn- 
chrotron spectrum 
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V{i>,<^f\\) = 



(B12) 



It is defined to be < 1 and ensure 77 — >■ 1 for 1/ <^ 
i/'c(i?max|| , Bo), as it is given by (|B11|) . In terms of i>, the 
modification factor is almost universal (i.e. allows for scal- 
ing with frequency). 

With the modification factor, the expression ()B8|) which 
describes the broad-band (radio-to- X-ray) synchrotron spec- 
trum from Sedov SNR becomes 



-(s-l)/2 



C(&,eK)r?(i>;ef||)ii:s||B^''+'''''i?'d-^(B13) 



F{u) = Cu 

The values of ^ are shown on Fig. IBll The parameter 
C is important in normalization of synchrotron spectrum: it 
varies in about 8 times over the parameter space. If injec- 
tion is considerably larger at parallel shocks (Ok < ^/3), 
the value of b is almost unimportant for amplitude of the 
synchrotron spectrum, but rather small changes in Ok may 
cause differences in in few times. In contrast, if injection 
tends to be isotropic (Ok > 2n/3), b plays the dominant 
role. 

In order to explore the parameter space, we made sev- 
eral runs to calculate the modification factors for different 
sets of parameters. Results are shown on Fig. IB2I where we 
also plot the experimental data in order to demonstrate rel- 
evance of the parameters for SN 1006. The modification fac- 
tor depends on ef y , b, Ok and s as well as on the function 
/e(0o). 

B2 Inverse-Compton emission 



The inverse-Compton flux (|B1|) from electrons in a black- 
body photon field with temperature Tcmb, at photon en- 
ergies far below TeV (i.e. when the Thomson regi me and 
power -law electron distribution are assumed, see IPetrukI 
l|2008l ) for details), is 



(B14) 



FT{e) = CTe~*'""/\T(6,OK)A:s||i?'d"' 
where e is the photon energy, 

2n IT 1 

Ct(6,Ok) = Jdip Jde sine fKjdaf^f.,K , (B15) 

00 

reflects the evolution of relativistic electrons downstream 
and 



os-l„2 ^ /I f„\(s+l)/2 /. \(s+5)/2 



where ctt is the Thomson cross-section, 



At(s) 



12 (s^ -I- 4s -I- 11) f z^-^+^y^dz 



t2 (s -I- 5)(s + 3)2 J exp{z) - 1 



2/(s + l) 



(B16) 



(B17) 



The contribution from electrons with energies around 
iJmax may be important for TeV 7-photons. The full expres- 
sion for IC process is 

Fic(e) = CicCic(e;6,OK,ef||,£n.ax||)A"s||ii'rf"' (BIS) 
where 

SaTkTcMBim^c^f-" 
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Figure B3. for different values of parameters h and 0k ■ s = 2 
1 




s, TeV 

Figure B4. Modification factor rji^. Lines are the same as on 
Fig. IB2I Experimental modification factor for SN 10 06 are shown 
for co mparison. It is obtained from the HESS data jAcero et al.l 
l20ld ) with the use of Eq. jlTJ and MF strength Bo = 3 /iG (upper 
panel) and Ba = 30 /^G (lower panel). 



2ft3c2 



(B19) 
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Cic(£; 6, Ok, ef||, -En 



= / d^p IdOsinOfK ldar^rs.K 



(B20) 



dll exp - 



7 



7n 



||f adf rad/E 



X(e,i5), 



where 7 is the electron Lorentz factor, T is an integral 
appearing in the theory of inverse-Compton process (e.g., 
IPetrukI boOSi ) : it accounts for the KN decline where neses- 
sary. 

In case s = 2 and Tcmb = 2.75, At = 0.710 and Ct = 
1.304 X lO^^^cgs, Cic = 1.186 x 10^^ cgs. 

In the limit E <^ -Emax and J? ^ et-Em ax, one has 
frad = 1 and X oc £, Bmin oc £^/^ l|Petrukll20ol ) and (fBlS)) 
transforms to (IB14[) . Therefore 

'(^-l)/2 



Cic(e) = CoE 



(B21) 



in this limit; Co — Ct/Cic. 

Let us introduce the modification factor for IC spec- 
trum: 



?)ic(e,ef||,Ejnax||) 



Cic(e, Efll , -Emax||)e 
CoCt 



(s-l)/2 



(B22) 



It is also defined to be rii^ < 1 and ensure r^ic — >■ 1 well 
below TeV energies. However, it is not so universal as for 
the synchrotron emission, Eq. (|B12|) : it does not scaled with 
the frequency and it depends on the absolute value of -Emax. 
The expression for the broadband IC spectrum is 

Eic(e) = CT£-^^-''/\T(6,eK)77ic(e,et||,i5^,,|l);f,l|i?^d-^(B23) 



The parameter behaves like ( (Fig. lB3 |l: it mostly de- 
pends on Ok for quasiparallel injection and on b for isotropic 
injection. However, the role of C^t is less important for nor- 
malization of IC spectrum because it varies in about 4 times 
over the parameter space. 

The modification factor of the IC specrum rjic is shown 
on Fig. IB4I in comparison with the observational data for 
SN 1006. It depends on efy, b, Ok, s and E^max as well as on 
the function /e(0o). 
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